Magic Wavelength for Atomic Motion Insensitive Optical Lattice Clocks 
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In a standing wave of light, a difference in spatial distributions of multipolar atom- field interactions 
may alter the definition of the magic wavelength to minimize the uncertainty of optical lattice clocks. 
We show that the magic wavelength can be determined so as to eliminate the spatial mismatch in 
electric dipole, magnetic dipole, and electric quadrupole interactions for specific combinations of 
standing waves. Experimental prospects of such lattices used with a blue magic wavelength are 
discussed. 
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Quantum absorbers trapped in well-designed electro- 
magnetic fields are the excellent candidates for future 
optical atomic clocks projecting uncertainties exceed- 
ing 10"-*^^ P, [2[, which represent the state of the art 
of the precision spectroscopy [1, 0]. While these atom 
traps provide long interrogation time and the Lamb- 
Dicke confinement of atoms necessary for ultrahigh res- 
olution spectroscopy, the relevant trapping fields impose 
an inherent fundamental limit on measurement uncer- 
tainties due to the atomic multipolar interactions with 
trapping fields [HI, 0, 0] and hyperpolarizability ef- 
fects J2|,l8|. It is of note that the detection and con- 
trol UIlO, [HI of the electric quadrupole interactions 
of ions with their trapping fields triggered an essential 
breakthrough for ion clocks operating on the S — D clock 
transitions to achieve the the uncertainties of 10~^^ and 
even below [1, [l2[ . 

Recently it was pointed out that the multipolar in- 
teractions of atoms with optical lattices may introduce 
spatial mismatch of the lattice potentials in the clock 
transition thus affecting optical lattice clocks' uncertain- 
ties at 10~^^ [13]. This inferred the slight breakdown of 
the original concept of the magic wavelength that cancels 
out the quadratic light shift in the clock transitions [2]. 
In this Letter, we discuss strategies to minimize the light 
shift uncertainties in optical lattice clocks by considering 
the electric dipole (El), magnetic dipole (Ml), and elec- 
tric quadrupole (E2) interactions of atoms in a standing 
wave. For specific lattice geometries, we show that the 
magic frequency ujm of a lattice clock can be defined so 
as to eliminate spatial mismatch of the lattice potentials, 
which is accompanied by a spatially constant differen- 
tial light shift Su. Consequently, optical lattice clocks 
free from atomic-motion-dependent clock shift are real- 
ized and the offset Su can be evaluated down to 10~^^ as 
a correction for the clock frequency. We discuss exper- 
imental prospects for the — ^Po clock transition of 
Sr atoms. In particular, combined with a blue-detuned 



magic wavelength [14], the proposed lattice geometry 
closely simulates the Paul trap employed in ion clocks 
[1], thus pushing lattice clocks' uncertainty towards the 
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regime. 



Optical lattices consist of a spatially periodic light 
shift formed by an interference pattern of electromagnetic 
waves. The simplest configuration is a one-dimensional 
(ID) lattice with a standing wave of light having wave 
vectors ±k and parallel electric fields E±k as shown in 
Fig. [TJ^a). The corresponding magnetic fields B±k are 
proportional to (±k) x E±k- As a result, anti- nodes 
of the total electric field E = E+k + E_k = 2E±k corre- 
spond to nodes of the magnetic field B = B+k+B_k = 0, 
and their amplitudes are a quarter of the wavelength 
A=27r/|k| out of phase, introducing different spatial de- 
pendences for the El and Ml interactions. Moreover, as 
the E2 interaction is proportional to the electric field gra- 
dient, its spatial dependence also differs from that of El. 
These different spatial dependences do not allow perfect 
cancellation of the quadratic light shift in the clock tran- 
sitions [13]. However, by admitting a constant differential 
light shift offset, we will show that spatial mismatch of 
the light shift can be eliminated, therefore the atomic 
motion dependent clock shift, which is detrimental to 
atomic clocks, can be removed. 

As there are hyperpolarizability effects [5, ^ that can- 
not be moderated by the magic wavelength, a blue de- 
tuned magic wavelength that confines atoms near the 
nodes of standing waves would be a promising choice, 
which will reduce the relevant uncertainties down to 
2 X 10"^^ [14]. Although the local electric field inten- 
sity for atoms is reduced, the E2 interaction, in turn, 
may severely affect the accuracy of a blue detuned lat- 
tice clock, as the electric field gradient can generally be 
maximum at the nodes of the standing wave. These blue 
detuned lattices can be realized using 3D lattice config- 
uration j^. Moreover, application of single-occupancy 
3D lattices with red and blue magic wavelengths may 
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effectively suppress the collision shift [16]. 

We consider 3D optical lattices consisting of three 
mutually orthogonal standing waves, whose results are 
straightforwardly applicable to lower dimensional lat- 
tices. In particular, we derive spatial dependences ^x(r) 
of the quadratic light shift with X = El, Ml, and E2 
interactions. For generality, we assume different electric 
field amplitudes = p^Eo of the three standing waves 
with equal wavelength A and wave number k=27r/X and 
different polarization vectors and (|p^| = |p^|=l) for 
the forward and backward running waves, respectively, in 
the ^=x, z directions denoted by unit vectors e^: 



E^(^, t) = E^ [p^ cos{k^ - cot) + p| cos{k^ + ut)] 
= E^ ^p^ cos cos ut + p^ sin k^ sin cot^ , 



(1) 



where we define p^=p^ ± p^. The total electric field 
vector in this lattice is given by E(r, t)=^^^^ y ^ E^(^, t). 

The principal contribution to the lattice potential is 
given by the second-order quasienergy shift due to the 
electric dipole atom-field interaction determined by the 
Hamiltonian VEi(r, Tg, t)=-d • E(r, [17], 



[/Ei(r) = -((V^|VEi5(re,t;r;,0VEi|^)) 
= — ^aEi(c^)(7Ei(r), 



(2) 



where Fg is the position vector of the outermost atomic 
electron relative to the atomic nucleus at r, d=— rg 
is the electric dipole moment, and ^(rg, t; r^, t') is the 
quasienergy Green function of an atom. Here the atomic 
units are used, e=m=h=l, where the speed of light is 
c^l37. The double angular brackets in Eq. ([2]) denote the 




FIG. 1: (a) Spatial distribution of an electromagnetic field for 
a ID standing wave, (b) Configuration of the electromagnetic 
fields and optical lattices for the cases (I)- (III), as described 
in the text. Optical lattice sites inside |x|, \y\, \z\ < 0.9A are 
indicated with their equipotential surfaces given by ^Ei(r) = 
0.3 and = 1. 



time integration (in variables t and t^) over the field os- 
cillation period T=2t:/uo and the spatial integration over 
the position Tg of the atomic electron. With the use of the 
electric field in Eq. ([1]), the spatial distribution function 
of the atom-field El interaction energy is given by 

^Ei (r) = - {^Yl P^^t + 2 ^^^^ y ' 

(3) 

To evaluate the contribution of the Ml interaction, it 
is sufficient to determine the magnetic field component of 
the lattice, which for each running wave with wave vector 
is given by B^(r, t)=k^/A:xE^(r, t). The total mag- 
netic field corresponding to the electric field in Eq. ([T]) is 
given by. 



B{r,t) 



E 

$,=x,y,z 



E^ (e^ X p^ sin /c^ sin cjt 

+ X p^ cos /c^ cos . (4) 

The magnetic dipole contribution to the lattice potential 
may be written similarly to Eq. ([2]), as the quasienergy 
shift corresponding to the atom-field Ml interaction is 
described by the Hamiltonian Vmi=— in • B(r, t), where 
m=— (J+S)/2c is the magnetic moment of an atom with 
atomic total momentum J and spin S. The spatial dis- 
tribution of the Ml interaction is given by. 



QMi{r) =i (^Y^ p^e^ X p+ sin k^^ 



\2 



(5) 



In the nonrelativistic approximation, the magnetic 
dipole polarizability in the n ^Sq ground state is zero, 
while for the n ^Pq met ast able state, it is given by 



<^Ml(^) = 



Er, 



3Pi 



■ Er, 



6c^[{E^sp^ - E^spJ^ 



(6) 



which is evidently the value of the second order in the 
fine-structure constant a=l/c. auii^) remains negative 
for CO higher than the n ^Pq — n ^Pi transition frequency. 

Not less important than the Ml Stark shift may be the 
contribution of the electric quadrupole (E2) interaction. 



f/E2(r) = -mVE2g{re,t;r'^,t')VE2m 
= — ^Q;E2('^)'2E2(r). 



(7) 



The value of the quadrupole polarizability aE2(^) is of 
the second order in the fine-structure constant a, just as 
the magnetic dipole polarizability in Eq. (|6|) is. The E2 
interaction operator may be taken from the Taylor series 
in powers of the small parameter /c|rg|<Cl for the total 
atom-electric-field interaction Hamiltonian, 



VE(r,rg,t) = i-g • ^ 



(re • V)^ 
(. + 1)! 



E(r,t), 



(8) 
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where all derivatives are taken with respect to the com- 
ponents of the position vector r, while is assumed to be 
constant. VE(r, re,t) includes all higher-order multipole 
interactions: The s=0 term corresponds to the Hamil- 
tonian VEi(r, and the 5=1 term to the Hamiltonian 
VE2(r, After substitution of this operator into Eq. (O 
and integration over time and angular variables, the spa- 
tial distribution of the quadrupole energy is determined: 

qE2{r) = \Y1 (^^^^ * + ^^^^ ' 

(9) 

where the sum runs over r]) = (x, 2:), and (2:, x). 

Correspondingly, the quadrupole polarizabilities of the 
ground and excited states are written in terms of the 
radial matrix elements, e.g., for the |?/^) = |n ^^o) state, 

"E2°(^) = ^{n'So\rligrD,+grnV^l\n'So), (10) 

where the radial Green functions of the singlet In- 
state subspace appear. 

Below we illustrate a few representative examples that 
allow cancellation of spatial mismatch of the lattice po- 
tentials. They are assorted by the electric field E/ and 
E5 of the forward and backward running waves that com- 
pose lattice standing waves, as summarized in Fig. [TJb). 

(I) E/IIE5 standing waves (p^=pp, in which we take 
Px=ey, Py=ez, and pz=^x- The El distribution is cal- 
culated to be, 

^Ei(r) = 2 (p^ cos^ kx + py cos^ ky + pi cos^ kz^ . (11) 

The Ml and E2 distributions are given by, 

<7Mi(r) = gE2(r) = Aq- gEi(r), (12) 

with Aq = 2{p1^p1^pI). Thus the distributions of Ml 
and E2 shifts in this lattice coincide and differ from g'Ei(r) 
in sign and by a constant offset of Aq. 

(II) Ej±E5 standing waves (p^±pp with their polar- 
ization vectors pointing at an angle 7r/4 to the standing- 
wave beams, i.e., Px = {ey^ez)/ V2, p^ = (-e^+e^)/y2, 
Py={ez^ex)/V2, p^y={-ez^ex)/V2, pz={ex^ey)/ V2 
and pl = {ex—ey)/^/2. The El and E2 distributions here 
coincide and Ml differs from them by sign and an offset, 

Qei (r) = qE2 (r) = Aq/2 + 2pxPz sin kx sin kz 

-\- 2pypz cos ky cos kz^ (13) 

(7Mi(r) = Aq- gEi(r). 

(III) E/±E5 standing waves (p^±p^), in which we take 
polarization vectors to be p^=e^±e2, p^=ez^ex^ and 



pf=ex^ey. Then the El and Ml distributions coincide 
and E2 differs from them by sign and an offset, 

<7Ei(r) =qMi{r) = Aq/2 ^ p^Py cos k{x ^ y) 

-^PyPz cos k{y + z) + pzPx cos k{z + x), (14) 
(7E2(r) =Aq - gEi(r). 

As indicated by these examples, it is essential that the 
spatial distributions of qx (r) may show the same spatial 
dependences apart from the sign and an offset Aq for 
particular lattice geometries. However, we note that this 
is not a general feature for optical lattices. For exam- 
ple, in the 3D lattice with E/IIE5 standing waves with 
Px=Pv=^z and Pz=^x employed in our previous experi- 
ment [16], neither g'Mi(r) nor qE2{^) shows the same spa- 
tial dependences as gEi(r), therefore the motional effects 
may limit clock uncertainties in future experiments. 

These three examples show prominent features by 
themselves. For a magic frequency, where the E2 (Ml) 
interaction is significantly larger than the Ml (E2) inter- 
action, case (II) [case (III)] will be more advantageous 
than the others, as the less significant Ml (E2) contribu- 
tion may well be neglected. For application to the blue 
magic wavelength, which highlights the reduction of the 
hyperpolarizability effects by trapping atoms near the 
nodes, case (I) would be a reasonable choice, as Eq. (pT]) 
suggests the creation of perfect nodes regardless of the 
intensity balance in the orthogonal lattice beams. Re- 
garding the lattice light polarization, case (I) shows lin- 
ear polarization, while cases (II) and (III) have elliptical 
polarizations that may give rise to the vector shift for 
atoms with non-zero angular momentum. 

Hereafter, we focus on case (I) and consider its appli- 
cation to the blue detuned magic wavelength. The clock 
transition frequency is expressed as 

i^ciock(^) = ^0 - ^AaEM(^)<7Ei(r)^o 

? (15) 
--Aao{io)AqEi^O{E^), 

where z/q is the atomic transition frequency, co is the lat- 
tice laser frequency. The quadratic light shift is decom- 
posed into two parts depending on their spatial depen- 
dences by AaEui^) = AaEi{uj) - Aauii^) ~ AaE2(^) 
and Aao{u;) = AaMi(^) + AaE2(^) using differential 
El, Ml, and E2 polarizabilities in the clock transition. 
Besides the hyperpolarizability effects that are minimized 
by use of a blue-detuned lattice [14], the magic frequency 
cOm may be given by AaEM(^m) = 0, allowing us to define 
LOrn independent of atomic motional states. The residual 
M1-E2 term Su = —^Aao{uJm)AqEQ provides an atomic- 
position-independent offset, which is solely related to the 
total lattice laser intensity AqE^ = 2{El^ El^ El). 

To find the magic frequency, it is essential to extract 
the atomic- motion-dependent term in Eq. (p!5|) . Harmon- 
ically approximating the trapping potential near the lat- 
tice node at x=y=z=\/ A and averaging over the atom 
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positions in the oscillator state |n) = \nx^ny^nz)^ the 
second term in the right-hand side of Eq. ([15]) should be 
replaced wth 



where 0.^^ =kp^EQ\J 2\a^l^{uj)\ / Ai is the vibrational fre- 
quency of atoms in the ^ direction of the lattice potential 
for the I = P(^Po) or S(^So) state. Here M is the mass 



of the atom, a^l^{u)^ 



«Ei(^)-«Mi(^)-^E2(^) is given 
by the El, Ml, and E2 polarizabilities in the £ state, and 
factor 2 accounts for the kinetic energy. 

Atomic-motion-dependent effects can be identi- 
fied by the clock frequency difference Au{u;^Sn) = 
z^ciock(^, (^n+n) — i^ciock(^, n) measured for atoms occupy- 
ing vibrational states differing by ^n. The magic fre- 
quency ujm can be determined by Ai/(cc;m, ^n)=0, 
equivalent to Aa^y^{uJm)=(^EMi^rn)-(^EMi^rn)=0. Once 
the magic frequency is determined, the residual M1-E2 
offset Su can be evaluated in terms of the vibrational fre- 
quencies of the lattice potential Q^=ft^ (^m)=^|^(^m) as 



6iy- 



MAao 
'2J? 



(17) 



Therefore all the essential measurements are done 
by the frequency measurements, once the magic fre- 
quency/wavelength is measured and shared. The same 
strategy should apply to ID optical lattices with red- 
detuned magic wavelength by setting px and py = 
= in case (I). 

The proposed optical lattice will be conveniently re- 
alized by a folded lattice jH, [l^, which maintains 
the relative phases of the orthogonal standing waves 
to realize linear lattice polarizations. Eq. ([T]) es- 
sentially assumes that the intensities of the counter- 
propagating beam pairs are balanced, which is accom- 
plished by preparing lattice beams inside an optical cav- 
ity |16|]. The blue-detuned magic wavelength for the 
Sr clock transition is experimentally determined to be 
Am=27rc/a;m^389.9 nm [14], where the numerical esti- 
mates for this lattice give Aao/|aEM|^— 1-4 • 10~^ with 
AaMi/AaE2^8 • 10"^ [20]. Therefore, the offset fre- 
quency is given by Su/27r^ 40-/ mHz for a trap frequency 
of ^}^/27T=75^/I kHz, where / is the single running wave 
laser intensity measured in kW/cm assuming p^ = 1. 
The uncertainty for this correction may be evaluated to 
be ~ 4/ mHz, assuming an inhomogeneity of the lattice 
intensity of 10 %. 

For ^^Sr atoms with a total angular momentum of 
F=9/2, the tensor light shift due to the spatial rotation 
of the lattice polarization with respect to the quantiza- 
tion axis may occur. As the shift is proportional to the 
light intensity near the nodes, the shift may be reduced 
to the mHz level for the blue magic wavelength. Bosonic 



isotopes such as ^^Sr or other atomic elements with nu- 
clear spin of 7=1/2, e.g., ^^""^Yb or ^^^Hg may well be 
used to suppress the tensor light shift. 

In summary, we present general formulae for the 
quadratic light shift taking multipolar atom-field inter- 
actions into account and show that the spatial mismatch 
of the interactions in the clock transition can be treated 
as a spatially constant offset Su for specific lattice ge- 
ometries. Numerical estimates are made for Sr atoms, 
and the relevant correction can be determined by trap 
frequency measurements with the mHz level. Combined 
with the blue magic wavelength, the hyperpolarizability 
effect is minimized, and clock uncertainty at the 10~^^ 
level will be within reach. 
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